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ABSTRACT 

The progenitors of Type la supernovae (SNe la) are not yet fully understood. The 
two leading progenitor scenarios are the single-degenerate (SD) scenario and the double¬ 
degenerate scenario. In the SD scenario, the collision of the SN la ejecta with its companion 
star is expected to produce detectable ultraviolet (UV) emission in the first few days after the 
SN explosion within certain viewing angles. A strong UV flash has recently been detected in 
an SN 2002es-like peculiar SN la iPTF14atg by Cao et al., which is interpreted as evidence of 
an early-time UV signature due to SN ejecta interacting with its companion star, supporting 
the SD scenario. In this paper, we present the expected luminosity distributions of early-time 
UV emission arising from SN la ejecta-companion interaction by performing binary popula¬ 
tion synthesis calculations for different progenitor systems in the SD scenario. Our theoretical 
predictions will be helpful for future early-time observations of SNe la to constrain their pos¬ 
sible progenitors. Assuming the observed strong UV pulse of iPTF14atg was indeed produced 
by the SN ejecta-companion interaction, our population synthesis model suggests that the pro¬ 
genitor system of iPTF14atg is most likely a red-giant donor binary system, and it is unlikely 
to have been a main-sequence or helium-star donor system. 
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1 INTRODUCTION 

Type la supemovae (SNe la) play a fundamental role in as¬ 
trophysics. The use of SNe la as cosmic distance indicators 
has led to the discovery of the accelerating expansion of the 
Universe (Riess et al. 1998; Perlmutter et al. 1999; Leibundgut 
2008). However, the nature of SN la progenitors and the 
physics of the explosion mechanism remain mysterious (see 
e.g. Hillebrandt & Niemeyer 2000; Maoz, Mannucci & Nelemans 
2014, for a review). It is generally accepted that SNe la arise from 
thermonuclear explosions of white dwarfs (WD) in binary sys¬ 
tems (Hoyle & Fowler 1960; Finzi & Wolf 1967; Nomoto 1982). 
Depending on the nature of the companion star, the most favored 
progenitor models of SNe la are classified into two general cat¬ 
egories, the single-degenerate (SD) scenario (e.g. Whelan & Iben 
1973) and the double-degenerate (DD) scenario (Iben & Tutukov 
1984; Webbink 1984). 

In the DD scenario, two carbon-oxygen (CO) WDs spiral in 
and merge due to gravitational wave radiation, resulting in a single 
object with a mass near the Chandrasekhar limit, which may then 
explode as an SN la (Iben & Tutukov 1984; Webbink 1984). Al¬ 
though the DD channel has been suggested to lead to an accretion- 
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induced collapse rather than a SN la (Nomoto & Iben 1985; 
Saio & Nomoto 1998), several hydrodynamical studies of mergers 
of DD systems have concluded that some DD pairs can explode 
as SNe la (Rasio & Shapiro 1995; Fryer et al. 2010; Pakmor et al. 
2010, 2011; Dan et al. 2011; Moll et al. 2014; Raskin et al. 2014). 
In particular, Pakmor et al. (2012) showed that the violent merger 
of two CO WDs with masses of 0.9 Mq and 1.1 Mq can lead to 
events that reproduce the observational characteristics of normal 
SNe la. 

In the SD scenario, a WD accretes material from its non¬ 
degenerate companion star, where the companion star could be 
either a hydrogen (H)-rich star such as a main-sequence (MS), 
a slightly evolved sub-giant (SG) star or a red giant (RG) star, 
or a helium (He) star. When the mass of the WD approaches 
the Chandrasekhar-mass limit, it explodes as an SN la (e.g. 
Whelan & Iben 1973; Nomoto 1982; Han & Podsiadlowski 2004). 
In the H-accretion scenario, only a fairly narrow range in the accre¬ 
tion rate above 10“^ Mq yr“' is allowed in order to attain stable H 
burning on the surface of the WD, avoiding a nova explosion. This 
makes the SD scenario difficult to explain the observed SN la rate 
(Han & Podsiadlowski 2004; Mannucci 2005; Wang et al. 2009; 
Wang & Han 2012; Maoz & Mannucci 2012). However, hydrody¬ 
namical simulations have shown that thermonuclear explosions of 
Chandrasekhar-mass WDs can reproduce some observational char- 
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acteristics of either subluminous Type lax supernovae (SNe lax) 
(Jordan et al. 2012b; Kromer et al. 2013, 2015; Fink et al. 2014), or 
normal SNe la (Gamezo, Khokhlov & Oran 2005; Seitenzahl et al. 
2013), or overluminous SNe la (Plewa, Calder & Lamb 2004; 
Meakin et al. 2009; Jordan et al. 2012a), which depend on the igni¬ 
tion conditions of a Chandrasekhar-mass WD. 

On the observational side, more and more evidence seems 
to favor the DD scenario, e.g., the non-detection of pre-explosion 
companion stars at the location of SNe la (Li et al. 2011; 
Bloom etal. 2012), the absence of H features in the nebu¬ 
lar spectra of SNe la (Leonard 2007; Lundqvist et al. 2013, 
2015; Shappee, Kochanek & Stanek 2013), the lack of radio 
and X-ray emission around peak brightness (Li et al. 2011; 
Bloom et al. 2012; Brown et al. 2012a; Chomiuketal. 2012; 
Horesh et al. 2012; Margutti et al. 2012), and the absence of a 
surviving companion star in SN la remnants (Kerzendorf et al. 
2009; Schaefer & Pagnotta 2012). All this evidence supports 
the DD scenario although most of them can also probably 
be explained by the so-called “spin-up/spin-down” SD model 
(Di Stefano, Voss & Claeys 2011; Justham 2011). In the “spin- 
up/spin-down” scenario, the fast spin of the WD can increase 
its critical explosion mass to be higher than the maximum mass 
achieved by the WD, leading to the WD must spin down before it 
can explode (Di Stefano, Voss & Claeys 2011; Justham 2011). As 
a result, the donor star might shrink rapidly before the WD explo¬ 
sion, because it would exhaust its H-rich envelope during a long 
spin-down timescale of the rapidly rotating WD until the SN la ex¬ 
plosion. In such cases, the companion star probably evolves to be 
a WD when the SN la explodes, which may explain the lack of H 
in late spectra of SNe la and the absence of a surviving compan¬ 
ion in the SN remnants (Di Stefano, Voss & Claeys 2011; Justham 
2011; Hachisu et al. 2012). On the contrary, some observations 
(Patat et al. 2007; Sternberg et al. 2011; Dilday et al. 2012) have 
detected signatures of H-rich circumstellar material (CSM) that are 
expected to exist around SNe la as the result of mass transfer from 
the companion, as well as WD winds (Nomoto 1982; Hachisu et al. 
1999). This supports the SD scenario. However, this detected abun¬ 
dant CSM is in conflict with the missing radio signal for other SN la 
events, e.g. SN 201 He (Chomiuketal. 2012; Horesh etal. 2012). 
Also, some studies suggest that CSM can also form in the DD sce¬ 
nario (Shen, Guillochon & Foley 2013; Soker et al. 2013). 

After the SN la explosion in the SD scenario, the SN ejecta 
expands freely for a few minutes to hours and then interacts with 
its non-degenerate companion star. As a result, the outer layers 
of the companion star are partially stripped and ablated while the 
star is shocked by the SN impact. Finally, the star survives the SN 
explosion and may show some peculiar features such as overlu¬ 
minosity or heavy element enrichment (Wheeler, Lecar & McKee 
1975; Marietta, Burrows & Fryxell 2000; Pakmor et al. 
2008; Liu et al. 2012, 2013a,b; Pan, Ricker & Taam 2012; 
Maeda, Kutsuna & Shigeyama 2014). The interaction of SN la 
ejecta with its companion star not only affects the star but also the 
SN itself. Kasen (2010) showed that the strong emission arising 
from the collision of SN la ejecta with its companion star in a 
SD progenitor system should be detectable in the first few hours 
to days after the SN explosion under favorable viewing angles. 
This strong emission is expected to alter the SN light curves at 
early times (Kasen 2010; Brown et al. 2012a; Cao et al. 2015; 
Moriya, Liu & Izzard 2015; Oiling et al. 2015). Compared to the 
brightness of the SN itself at early times, such strong emissions 
should be brightest in the ultraviolet (UV). Therefore, early UV 
observations are proposed as a direct way to test progenitor models 



Figure 1. Distributions of binary separations at the moment of SN explosion 
in different SD progenitor scenarios. 


of SNe la. Following the analytical method of Kasen (2010), 
analysis of observed early light-curves of SNe la has been carried 
out by different groups to look for evidence of shock emission 
and thus constrain possible progenitors of SNe la (Hayden et al. 
2010; Brown et al. 2012a,b; Cao et al. 2015; Oiling et al. 2015). 
Most interestingly, a strong UV flash after the SN explosion has 
recently been detected from early-time observations of a peculiar 
subluminous SN la iPTF14atg by Cao et al. (2015). They further 
interpreted this early UV flash as strong evidence of excess 
luminosity produced by SN ejecta interaction with its companion 
star, and suggested that some SNe la arise from the SD scenario 
(Cao et al. 2015). 

In this work, we use binary population synthesis (BPS) mod¬ 
els for different progenitor systems of SNe la in the SD scenario 
to obtain the progenitor properties at the moment of SN la ex¬ 
plosion, e.g., the companion radius (R) and the binary separation 
(a). We then put these progenitor parameters into the analytical 
model of Kasen (2010) to predict the distributions of early UV lu¬ 
minosity arising from the SN ejecta interacting with its compan¬ 
ion star. Comparing our theoretical predictions with future early 
observations of SNe la will be helpful to constrain their possible 
progenitors. In addition, assuming the strong UV flash detected in 
iPTF14atg was indeed produced from the SN ejecta-companion in¬ 
teraction, we put further constraints on the progenitor system of 
iPTF14atg based on our BPS results. 

The paper is organized as follows. In Section 2, we describe 
the BPS method used in this work. Theoretical early UV emission 
of SN ejecta-companion interaction are calculated in Section 3. In 
Section 4, we compare our predictions with the early-time observa¬ 
tions of iPTF14atg and some discussions are presented. Our con¬ 
clusions are summarized in Section 5. 


2 BINARY POPULATION SYNT HE SIS CALCULATIONS 

Adopting the method described in Han & Podsiadlowski (2004), 
we have performed BPS calculations for the SD Chandrasekhar- 
mass scenario of SNe lax by considering different types of com¬ 
panion stars in Liu et al. (2015). With that BPS calculation, we 
have obtained the distributions of progenitor properties of differ¬ 
ent SD binary systems at the moment of SN explosion although we 
assumed all Chandrasekhar-mass CO WDs would lead to weak de¬ 
flagration explosions and hence to SNe lax (Kromeretal. 2013; 
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Figure 2. Distributions of theoretical light curves predicted from the interaction between SN la ejecta and its MS companion star in the four UV bands (wvwl, 
uvm2, uvm2 and U) of the Swijt telescope. Here, 6 is the viewing angle, the maximum UV flux corresponds to a viewing angle of 0 = 0° (looking down on the 
companion star). For a given viewing angle, each shaped region shows the covered magnitude range of our theoretical predictions with the MS donor scenario 
at different epochs after the SN explosion. We assume the ejecta mass Mej = 1.4 Mq and the SN explosion energy £ej = 1.0 X 10^* erg. 


Finketal. 2014) in that work. In the SD scenario, whether a 
Chandrasekhar-mass CO WD would finally lead to a sublumi- 
nous (i.e., lax events), normal, or overluminous SN explosion, 
depends on the exact ignition conditions of the CO WD. How¬ 
ever, as the WDs in our previous one-dimensional binary evolu¬ 
tion calculations (Liu et al. 2015) were treated as point masses, 
it is impossible to determine the exact ignition conditions of the 
Chandrasekhar-mass CO WD. Therefore, we cannot distinguish 
whether the Chandrasekhar-mass CO WD undergoes a delayed- 
detonation explosion that matches normal SNe la or a weak de¬ 
flagration explosion of SNe lax, and the progenitor properties at 
the moment of SN explosions are the same in our BPS calcula¬ 
tions for the SD Chandrasekhar-mass scenario. Therefore, we di¬ 
rectly use the data obtained from BPS calculations of Liu et al. 
(2015) to calculate early UV emission of shocked ejecta of SNe 
la in the SD scenario in this work. Here, we only briefly de¬ 
scribe our BPS method and basic assumptions, and we refer to 
Han & Podsiadlowski (2004) for a detailed description about the 
BPS method used in this work. 

We use the Cambridge stellar evolution code STARS 
(Eggleton 1971, 1972, 1973) to trace the detailed binary evolution 
of a set of systems consisting of a CO WD and a MS or SG (which 
is called the MS donor scenario in the following), a RG, or a He 
star companion. Once the companion star fills its Roche lobe, the 
mass transfer occurs though Roche-lobe overflow (RLOF) which 
is treated by the prescription of Han, Tout & Eggleton (2000). 
Once the WD has accreted the transferred H-rich or He-rich ma¬ 


terial from the companion star and increased its mass to near 
the Chandrasekhar-mass limit (~ 1.4 M©), we then assume an 
SN la explosion occurs. The optically thick wind assumption 
of Hachisu, Kato & Nomoto (1996) is used to describe the mass 
growth of a CO WD by the accretion of H-rich material from the 
donor star. The prescription of Kato & Hachisu (2004) is imple¬ 
mented for the mass accumulation efficiency onto the CO WDs 
when He-shell flashes occur. With a series of binary evolution cal¬ 
culations for various close WD binary systems, we determine the 
initial parameter space leading to SNe la in the orbital period¬ 
secondary mass (i.e. logio ) plane for various initial CO WD 

masses (Af^^p) in the MS, RG, and He star donor channel. 

Eurthermore, to obtain the distributions of progenitor prop¬ 
erties of binary systems at the moment of SN explosion in 
the SD scenario, we use a rapid population synthesis code 
(Hurley, Pols & Tout 2000; Hurley, Tout & Pols 2002) to perform 
a detailed Monte Carlo simulation. When a binary system evolves 
to a CO WD 4- MS (SG, RG, or He star) system and achieves 
the beginning of the RLOF phase, whether this binary system can 
lead to a SN la explosion is determined by checking if it is lo¬ 
cated in the SN la production regions in the plane of (logm P', 
M'p for its based on our detailed binary evolution. Then, 

progenitor properties at the moment of SN explosion are obtained 
by interpolation in the three-dimensional grid (M(yp, M^, logic P') 
of the close WD binaries obtained in our detailed binary evolu¬ 
tion calculations. In our BPS calculations, the initial mass func¬ 
tion of Miller & Scalo (1979) is used. We assume a circular bi- 
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Figure 3. As Fig. 2, but for the RG donor Chandrasekhar-mass scenario. 


nary orbit and set up a constant initial mass ratio distribution (i.e. 
n(q') =constant, see Goldberg & Mazeh 1994; Bender & Simon 
2008; Duchene & Kraus 2013). The initial separation distribution 
of Han, Podsiadlowski & Eggleton (1995) is adopted. We assume a 
constant star formation rate (Han 2008). The standard energy equa¬ 
tions of Webbink (1984) are used to calculate the output of the CE 
phase. The CE ejection is determined with two highly uncertain pa¬ 
rameters, acE and T. Here, acE is the CE ejection efficiency, i.e. the 
fraction of the released orbital energy used to eject the CE and A is 
a structure parameter that depends on the evolutionary stage of the 
donor star. In this work, the models with a metallicity of Z = 0.02 
and a parameter of acET^O.S are adopted (Han & Podsiadlowski 
2004). 


3 PREDICTIONS OF THE EARLY UV EMISSION 

In the analytical model of Kasen (2010), SN ejecta collides with 
its companion star, the impacting layers are re-shocked and the 
kinetic energy partially is dissipated, causing UV flux emission 
from the shock-heated region. This emission dominates the early- 
time (few days after SN explosion) light curves of SNe la within 
certain viewing angles (Kasen 2010). At later epochs, this shock 
flux becomes weak and the light curves of SNe la are dominated 
by the flux of the SN itself. Kasen (2010) analytically estimated 
the early isotropic bolometric luminosity Lc.iso by the collision as 


(Equation 22 of Kasen 2010) 



where a is the binary separation at the moment of SN explosion, 
Mej is the SN ejecta mass, Vej is the velocity of the SN ejecta collid¬ 
ing with the binary star, aTc is the electron scattering opacity in the 
SN ejecta which is assumed to he a constant value in this work, and 
1 is the time since SN explosion. This luminosity is strongly depen¬ 
dent on the binary separation at the moment of SN la explosion. In 
addition, the effective temperature of the prompt shock emission is 
estimated as follows (Equation 25 of Kasen 2010). 


4 

The strong UV emission at early times from the SN ejecta in¬ 
teracting with the companion star not only strongly depends on the 
progenitor system but also on the viewing angle 6 (Kasen 2010). 
Eor a given binary progenitor system, the maximum excess emis¬ 
sion should be observed when the viewing angle is 0°, i.e., the 
companion star lies directly along the line of sight between the ob¬ 
server and the SN explosion. On the contrary, the excess emission 
detected for viewing angle 6 = 180° should be negligible, corre¬ 
sponding to a geometry in which the SN la lies directly in the line 
of sight between the observer and the companion star. Brown et al. 
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Time since expiosion (days) Time since expiosion (days) 

Figure 4. As Fig. 2, but for the He star donor Chandrasekhar-mass scenario. 


(2012a) estimated the dependence of the early emission from SN- 
companion interaction on the viewing angle. They found that the 
dependence of fractional flux values (/) on the viewing angle (6) 
can be fitted by a function (Equation 3 of Brown et al. 2012a) of 
the form: 

/ = (0.5 cos B + 0.5) X (0.14 0^ - 0.4 0+1), (3) 

where the viewing angle 0 is in unit of radians. In this work, we di¬ 
rectly adopt this function to calculate the observed flux at different 
viewing angles for a given epoch. 

To facilitate direct comparison with the early-time obser¬ 
vations of SNe la, the early isotropic bolometric luminosity of 
shocked gas is converted to broad band magnitudes by using the 
same method as Brown et al. (2012a). Adopting the progenitor 
properties at the moment of SN la explosion obtained from our 
BPS calculations in Section 2, we create a set of blackbody mod¬ 
els under a given viewing angle with appropriate temperature and 
luminosity based on Equations 1, 2 and 3. In Fig. 1, we present the 
distributions of binary separation at the moment of SN explosion in 
our BPS calculations for three different SD progenitor scenarios. It 
is shown that most progenitor systems in the He star donor scenario 
have a relatively small separation (less than 10** cm) at the moment 
of SN explosion, and the binary systems in the H-rich donor sce¬ 
nario cover a separation range of about 10** - 4.0 x 10*^ cm. With 
the created blackbody models, we then calculate the synthetic mag¬ 
nitudes at different bands by using transmission curves of chosen 
filters. Specifically for this work, the 5wi/T/UVOT filter system and 
their corresponding transmission curves are selected. The absolute 
magnitudes are calculated in the AB magnitude system. 

With different donor star models, the theoretical distributions 


of UV emission at early times from SN ejecta-companion interac¬ 
tion in the SD scenario are presented in Figs. 2, 3 and 4. Here, 
we assume typical explosion parameters of an SN la, i.e., the 
amount of ejecta mass = 1.4 A/q and the SN explosion energy 
£ej = 10 X 10^* erg. These typical values lead to a mean expansion 
velocity of SN ejecta of 10"'kms“*. Here, only four different 
SwiftfnWOT filters are chosen to show as examples, i.e., the uvwl, 
uvm2, uvm2 and U bands (the B and V band are not included), be¬ 
cause the excess emission at early times should be brightest in the 
UV and become subordinate at longer optical wavelengths (Kasen 
2010; Brown et al. 2012a). For a given viewing angle, we show the 
upper-limit and lower-limit absolute magnitude in each band as a 
function of time since the SN explosion. It is shown that the the¬ 
oretical UV luminosities of shocked gas in the MS donor scenario 
cover a relatively wider magnitude range because binary systems 
in this scenario have a wide range of progenitor properties at the 
moment of SN explosion compared to those in other two donor sce¬ 
narios (Fig. 1 ). In addition, the early UV luminosities in the He star 
donor scenario are much weaker than those of the H-rich scenario 
by more than three magnitudes because He star donor progenitor 
systems have closer separations (Fig. 1) when the WDs grow to 
the Chandrasekhar-mass limit (Fig. 4). We note that the explosion 
parameters of typical SNe la used here may not match the case of 
iPTF14atg because it is identified as a subluminous SN 2002es- 
like event (see Section 4.1), which however, strongly depends on 
the exact explosion mechanism of a Chandrasekhar-mass CO WD. 

In Figs. 2-4, we present possible magnitude ranges that can 
be covered by the early UV luminosities arising from the ejecta- 
companion interaction in different SD progenitor scenarios. To bet¬ 
ter compare our predictions with early-time UV observations of 
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Figure 5. Left: the probability distribution of the absolute magnitudes of shocked gas in uvm2 band at different epochs after the SN la explosion for the MS 
(top), RG (middle) and He star (bottom) donor scenario, respectively. The fitting curves of the distributions are also plotted except for the RG donor scenario. 
Right: comparison between the early observations of a peculiar SN la iPTF14atg and our predictions. The boundaries (short thin dashed curves) of gray region 
con'espond to the upper-limit and lower-limit magnitude in uvm2 band. The long bold dashed lines show absolute magnitudes at the peaks of fitting curves of 
corresponding figures in left panel. The location of the early UV flash detected in iPTF14atg (Cao et al. 2015) is shown with the square in red simply assuming 
10% error of the time because of the uncertainty of the explosion date. 


SNe la and put stronger constraints on the possible progenitors 
of SNe la, the probability distributions of our predicted UV lumi- 
noisties at a given epoch after the SN explosion are further shown 
in Fig. 5. Here, we only consider the case with a viewing angle 
9 = 0°, i.e., the maximum flux case, and the results at uvml band 
are chosen to give an example because uvwl and uvw2 filters have 


an extended red tail. The distributions of different viewing angles 
and other bands can be roughly estimated according to the results in 
Fig. 2 to Fig. 5. Also, only the results within the first five days after 
the SN explosion are presented in Fig. 5 because SN light curves 
at late times are dominated by the flux from the SN itself (Kasen 
2010). As it is shown in Fig. 5, a peak is clearly seen in the proba- 
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bility distributions of the MS and He star donor scenario. However, 
the probability distributions of absolute magnitudes within a given 
day in the RG scenario is relatively flat (Fig. 5). This may indicate 
that early UV luminosities of different progenitor systems in the 
RG donor scenario have a similar probability at a given epoch after 
SN explosion. 

4 DISCUSSIONS 

4.1 Comparison with the early UV flash of iPTFI4atg 

A strong UV pulse in the early-time light curves of a subluminous 
SN la iPTF14atg has been detected by Cao et al. (2015). They fur¬ 
ther suggested that this strong UV pulse supports that some SNe la 
arise from the SD Chandrasekhar-mass scenario. Comparing their 
observations to analytical models of Kasen (2010), they found that 
the theoretical UV emission from shocked ejecta of a binary sys¬ 
tem with a separation of about 60 Rq can match the observations 
assuming the explosion energy £ej = 10^* erg and ejecta mass 
Mej = 1.4 Mq. Here, we also compare the predictions from our 
BPS calculations for different SD scenarios to the early UV flash 
of iPTF14atg to put a more stringent on the possible progenitor of 
this event. 

A detailed comparison with the results of the uvm2 band is 
given in the right column of Fig. 5. It shows that the observed early 
UV pulse of iPTF14atg is stronger than the upper limit of the the¬ 
oretical UV luminosity predicted from the MS donor and He star 
donor scenario at the same epoch (Fig. 5). Considering that the ac¬ 
tual viewing angle may be larger than 6 = 0° and taking the uncer¬ 
tainties on the explosion parameters (see Section 4.4) into account, 
it is unlikely that the detected UV flash in early-time light curves of 
iPTF14atg arises from the interaction of SN la ejecta with a MS or 
He star companion. However, our theoretical early UV luminosi¬ 
ties from the RG donor scenario can include the early UV flash 
of iPTF14atg (Fig. 5) although the theoretical upper-limit luminos¬ 
ity in this scenario can still decrease to a value weaker than the 
observed UV flash in iPTF14atg as the viewing angle increases. 
Therefore, whether the predictions from the RG donor scenario can 
support the early UV observations of iPTF14atg depends on the ac¬ 
tual viewing angle of the observations. Nevertheless, we still sug¬ 
gest that the RG donor binary system as progenitor of iPTF14atg is 
more possible if the early UV flash detected in an SN 2002es-like 
SN iPTF14atg indeed arises from the ejecta-companion interaction 
in the SD scenario. 


4.1.1 Different explosion parameters 

The comparisons given above are based on an assumption of the 
typical SN la explosion parameters as those used in Cao et al. 
(2015). However, it has been suggested that iPTF14atg is more 
likely to belong to the subluminous SN 2002es family by Cao et al. 
(2015). Meanwhile, it is shown that the peak magnitude of 
iPTF14atg is similar to SN 2005hk (a typical SN lax event). The 
hydrodynamical simulations of Kromeretal. (2013) have shown 
that the off-centre-ignited weak deflagration of Chandrasekhar- 
mass CO WDs is able to reproduce the characteristic observational 
features of SN 2005hk quite well. Adopting the weak deflagration 
explosion model (£ej = 1.34 x 10^° erg and Mej = 0.37 Mq, which 
corresponds to a mean expansion velocity Vej ~ 6 000 km s“') of 
Kromeretal. (2013) as inputs of the analytical model of Kasen 
(2010), we predict theoretical early UV luminosities in the uvml 



Figure 6. Similai' as Fig. 2, but for theoretical predictions of SN lax case. 
Here, only a viewing angle 8 = 0° and uvml band are considered. The ejecta 
mass (Mej = 0.37 Mq) and SN explosion energy (£ej = 1.34 X 10^® erg) are 
directly derived from a weak deflagration explosion of the Chandrasekhar- 
mass C/0 WD for SN lax event (Kromer et al. 2013). The location of early 
UV flash of iPTF14atg is also shown (Cao et al. 2015). 


band of different SD progenitor scenarios (Fig. 6). Compared to the 
results with typical explosion parameters in Figs. 2-4, it is shown 
that the early UV luminosities of shocked ejecta are weaker by 
about one magnitude when the explosion parameters of lax events 
are used (Fig. 6). If the SN iPTF14atg was produced from an ex¬ 
plosion model that has similar explosion parameters with the weak 
deflagration explosion model for SNe lax in Kromer et al. (2013), 
the observed UV flash in iPTF14atg is a bit stronger than the upper- 
limit luminosity of our RG donor model at similar epoch (Fig. 6). 

However, the actual explosion date was not well constrained 
in Cao et al. (2015). Therefore, the delay time between the SN ex¬ 
plosion and the moment of ejecta-companion interaction happens 
is not well determined, which means that the location of early UV 
pulse of iPTF14atg in Figs. 5 and 6 can be shifted a bit towards 
the left (or the right) if the actual explosion date is earlier (or later) 
than that used in Cao et al. (2015). 

4.1.2 SN birth rate 

In our BPS calculations, we obtain a SN la rate of about 3 x 
10“^ yr“' (i.e., three SN la events every 10^ years) in the RG donor 
channel for a constant star formation rate (5.0 MQyr^*), which 
is about 1% of the total SN la rate if we assume that the total 
SN la rate is 3.5 x 10“^ yr“' based on the inferred Galactic SN 
la rate of 3 - 4 x 10“^ yr“' (Cappellaro & Turatto 1997). Obser- 
vationally, Ganeshalingam et al. (2012) estimate that roughly 2.5% 
of SNe la may belong to SN 2002es-like events within a fixed vol¬ 
ume. The birth rate from the RG donor channel in our BPS calcula¬ 
tion is roughly comparable with observed rate of SN 2002es-like 
event (but see Hachisuetal. 1999). In addition, the delay times 
between star formation and SN explosion in the RG donor sce¬ 
nario are quite long (> 3 Gyr), which seems to support the fact 
that 2002es-like SNe la are likely from an old stellar population 
(Ganeshalingam et al. 2012) and the host galaxy of iPTF14atg is 
classified as an old galaxy. 

However, the contribution of the RG donor scenario to 
the total SN la rate is quite uncertain (Hachisu et al. 1999; 
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Figure 7. Simitar to Fig. 5, but showing the results with a random distribution of viewing angle 6. The distributions by adopting the explosion parameters for 
typical SNe la (left column) and SN lax event (right column) are presented. 


Han & Podsiadlowski 2004; Chen, Han & Tout 2011). If mass 
transfer in symbiotic binaries (RG donor model) occurs through 
the stellar wind rather than the RLOF assumed in this work 
(Branch et al. 1995), the SN contribution from this scenario could 
remarkably increase and binary systems could have much larger 
orbital separations than those in this work at the moment of SN la 
explosion (Chen, Han & Tout 2011). Subsequently, early UV emis¬ 
sions from shocked gas may be much stronger than those obtained 
from our RG donor models. 


4.2 The elTect of viewing angle 

In Fig. 5, we present the expected magnitudes in the uvml band 
from a fixed viewing angle (0 = 0°). However, we do not neces¬ 
sarily observe the interaction from 0 = 0°. To investigate the ef¬ 
fect of the viewing angle on the observational magnitude, we ran¬ 
domly assign the viewing angle and obtain the probability distri¬ 
bution for the observational magnitude in the uvml band. Fig. 7 
shows the probability distribution for both normal SNe la (i.e. 
£ej = 1.0 X 10^* erg and Mej = 1.4 M©) and SN lax events (i.e. 
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£ej = 1.34 X 10^® erg and - 0.37 Mq). Comparing to the dis¬ 
tributions in Fig. 5, the peak of the distribution becomes about one 
magnitude fainter after taking the effect of the viewing angle into 
account. However, the maximum magnitude achieved by the SN 
ejecta-companion interaction does not change and our conclusion 
regarding the SN iPTF14atg remains unaltered. 

4.3 The Spin-up/Spin-down model 

In the SD scenario, a WD accretes and retains companion mat¬ 
ter that carries angular momentum. As a consequence the WD 
spins with a short period which leads to an increase of the crit¬ 
ical explosion mass. If the critical mass is higher than the ac¬ 
tual mass of the WD, the SN explosion can only occur after 
the WD spins down to decrease its critical explosion mass to 
fall below the current mass of the WD with a specific spin- 
down fimescale (Di Stefano, Voss & Claeys 2011; Justham 2011; 
Hachisu et al. 2012). In such case, the MS or RG donor star might 
shrink rapidly before the SN la explosion occurs because most its 
H-rich envelope has been exhausted during a long spin-down phase 
of the rapidly rotating WD. As a result, the early UV signature from 
SN la ejecta interacting with its companion star should be remark¬ 
ably reduced compared to the results we presented in this work 
because the donor star is much smaller than its Roche lobe at the 
moment of SN la explosion. In addition, as previously mentioned, 
this “spin-up/spin-down” model seems to be able to explain tbe lack 
of H in late spectra of SNe la and possibly the absence of a surviv¬ 
ing companion star in the SN remnants (Di Stefano, Voss & Claeys 
2011 ; Justham 2011 ; Hachisu et al. 2012). However, the exact spin- 
down timescale of the WD in this model is quite unknown. 

4.4 The effect of different explosion models 

In the SD scenario, depending on the exact ignition con¬ 
dition, a Chandrasekhar-mass CO WD can undergo a de¬ 
flagration, a detonation, or a delayed detonation explosion 
(Arnett 1969; Woosley, Taam & Weaver 1986; Khokhlov 1989; 
Rbpke & Hillebrandt 2005; Rbpke et al. 2007), which could finally 
produce SN lax events (Jordan et al. 2012a; Kromeretal. 2013, 
2015; Fink et al. 2014), normal SNe la (Seitenzahl et al. 2013), or 
overluminous SNe la (Plewa, Calder & Lamb 2004; Jordan et al. 
2012b). The distributions of early UV emission of shocked gas in 
normal SNe la and lax events have been presented in Sections 3 
and 4.1 . 

By analytically deriving the existence of a characteristic 
length scale which establishes a transition from central ignitions 
to buoyancy-driven ignitions, Fisher & Jumper (2015) suggested 
that the Chandrasekhar-mass WDs in the SD scenario are generally 
expected to explode as overluminous SNe la rather than normal 
SNe la or lax events. If the suggestion of Fisher & Jumper (2015) 
is true, early UV signatures due to the ejecta-companion interac¬ 
tion may be easier to be detected in overluminous SNe la because 
overluminous SNe la may have a higher explosion energy than that 
of normal SNe la. A higher explosion energy is expected to cause 
stronger early UV emissions of shocked gas for a given binary pro¬ 
genitor. 

4.5 Model uncertainties 

The prescription of the optically thick wind model of 
Hachisu, Kato & Nomoto (1996) and He-retention efficien¬ 


cies from Kato & Hachisu (2004) are used to describe the mass 
accumulation efficiency of accreting WDs in our detailed bi¬ 
nary evolution in this work. However, different mass-retention 
efficiencies may lead to the results that are different from our 
BPS calculations (e.g., see Bours, Toonen & Nelemans 2013; 
Piersanti, Tornambe & Yungelson 2014; Ruiteretal. 2014; 
Toonen et al. 2014). However, the exact mass-retention efficiency 
is still not well-constrained. 

It should be kept in mind that the initial conditions of BPS 
calculations such as the CE evolution, current star-formation rate 
and initial mass function, may be sensitive to the assumed parame¬ 
ters in specific BPS codes, which lead to some uncertainties on the 
BPS results (Toonen et al. 2014). However, current constraints on 
these parameters (e.g., the CE efficiency, see Zorotovic et al. 2010; 
De Marco et al. 2011 ; Ivanova et al. 2013) are still weak. Eor a de¬ 
tailed discussion for the effect of different theoretical uncertainties, 
see Claeys et al. (2014). 

Einally, we point out that the analytical model of Kasen (2010) 
is simple, which may lead to the predicted early UV luminosity 
being overestimated, or underestimated. Eor instance, a constant 
opacity from electron scattering is assumed in the analytical model 
of Kasen (2010), the luminosity can be reduced if the realistic opac¬ 
ity in the interacting material is decreased (Brown et al. 2012a). To 
put stronger constraints on the possible progenitor of SNe la, the ra¬ 
diation hydrodynamical simulations of the interaction between SN 
la ejecta and its companion star are encouraged to address in the 
future study. 


5 CONCLUSION AND SUMMARY 

In the SD Chandrasekhar-mass scenario of SNe la, strong UV emis¬ 
sions arising from the SN ejecta interacting with its companion 
star is suggested to be dominant in the early-time light curves of 
SNe la (Kasen 2010). These prompt UV emission characteristics 
are expected to constrain possible progenitor systems of SNe la. 
In this work, we perform BPS calculations for different progeni¬ 
tors of SNe la within the SD Chandrasekhar-mass scenario to ob¬ 
tain progenitor properties at the moment of SN la explosion. We 
then predict the distributions of early-time UV luminosities of SN 
ejecta-companion interaction of different progenitors by combing 
our BPS results with the analytical model of Kasen (2010). These 
theoretical predictions for early UV emissions of shocked ejecta 
will be helpful to compare with future early-time observations of 
SNe la to put constraints on their possible progenitors. 

In addition, assuming the observed strong UV pulse of 
iPTF14atg indeed arises from the SN ejecta-companion interaction, 
we have compared our theoretical UV luminosities to the observed 
UV pulse at early times of a peculiar subluminous SN la iPTP14atg 
(Cao et al. 2015). With our BPS models, we find that predicted UV 
luminoisties in the MS and He star donor scenario struggle to reach 
an upper-limit absolute magnitude as strong as early UV flash ob¬ 
served in iPTP14atg, which suggesfs fhaf the MS and He star donor 
binary systems are unlikely to be the progenitor. In particular, the 
He star donor binary system as progenitor of SN iPTP14atg can be 
ruled out by our population synthesis calculations. However, the¬ 
oretical UV luminosities at early times in our RG donor scenario 
provide a match with the early UV observations of iPTF14atg un¬ 
der some appropriate conditions such as a favorable viewing angle. 
Therefore, we suggest that a red-giant donor binary system as the 
progenitor of an SN 2002es-like peculiar SN la iPTP14atg is more 
possible based on our populations synthesis model. 
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